Cosmological viability of f(T) models according to Noether symmetries 
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Noether's theory offers us a useful tool to research the conserved quantities and symmetries 
of the modified gravity theories, among which the f(T) theory, a generally modified teleparallel 
gravity, has been proposed as an alternative gravity model to account for the dark energy 
phenomena. By the Noether symmetry approach, we investigate the power-law, exponential and 
polynomial forms of f(T) theories with corresponding conserved quantities and symmetries. All 
forms of f(T) concerned in this work possess the time translational symmetry, which is related 
with energy condition or Hamilton constraint. And we find out that the performances of the 
power-law and exponential forms are not pleasing. It is rational adding a linear term T to T n as 
the most efficient amendment to resemble the teleparallel gravity or General Relativity on small 
scales, ie., the scale of the solar system. The corresponding Noether symmetry indicates that 
^»i ■ only time translational symmetry remains. Through numerically calculations and observational 

*—* • data-sets constraining, the optimal form aT + fiT^ 1 is obtained, whose cosmological solution 

resembles the standard ACDM best with lightly reduced cosmic age which can be alleviated by 
introducing another T m term. Proper constraints on the ratio /3/ct are also given, which can 
be converted into constraints on the ratio of the magnitude of torsion between large and small scales. 
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I. INTRODUCTION 



> 

00 ' Diverse sets of astrophysical observations, which include WAMP 1], SN la 0, X-ray Q, LSS 0] and SDSS [| 

lO ! clearly indicate that our universe is currently accelerated expanding. These observations also imply that our universe 
^*>0 is spatially flat on large scales, and consists of about 68% dark energy (DE) with effectively negative pressure, 32% 
tH- dust matter (cold dark matter plus baryons), and negligible relativistic constituent (including radiation and neutrinos) 

f^*) ■ at present. 

In order to explain the cosmic acceleration, many theories have been proposed phenomenologically and physically. 
Generally speaking, approaches can be classified into two broad groups: Modification the energy-momentum tensor 
J> ' of Einstein's equation, where the framework of general relativity (GR) is kept unchanged and specific dark energy 
fields are invoked. The simplest candidate of dark energy is a tiny positive time-independent cosmological constant 
rS ■ A with equation-of-state (EOS) parameter u=— 1. Although this is highly consistent with the available observational 
data-sets, it faces difficulties such as the microphysical origin, besides, recent observations favor EOS parameter of 
constant DE lower than -1 [6j. Modified gravity theories provide another way, which present natural unification of 
the early-stage in action and late-stage acceleration due to different roles of gravitational terms relevant at small 
and at large curvature. Scalar-tensor theories [7J], f(R) gravity [SHTlj. and others are studied extensively. Besides 
which the f(T) (l2T - tl4j gravity, where T is the so-called torsion scalar, has been proposed recently and attracted 
much attention. In the case f(T)=T, f(T) theory can be directly reduced to the Teleparallel Equivalent of General 
Relativity (TEGR) [15|,[l6| which was first propounded by Einstein in 1928 [17| to unify gravity and electromagnetism 
by introducing a vierbein field along with the concept of absolute parallelism or teleparallelism. 

Constraints on f(T) gravity [18| by latest observational data-sets, dynamical behavior [lj| and cosmic large scale 
structure [20] , relativistic neutron star [21] , matter bounce [22[ and perturbations 23] in f(T) gravity framework, 
static spherical symmetry solutions, validity of Birkhoff 's theorem [241 ] and the EOS parameter crossing the phantom 
divide [251 ] are investigated by previous researches. 
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Without fundamental principles, most of the models generated from modified gravity theories in previous work 
were proposed phenomcnologically, which seem arbitrary. Specific solutions of the field equations depend on explicit 
expressions of MG (modified gravity) theories, which are important for gaining insight into the mathematical and 
physical contents of the theories. In our present work, we concern about whether there exist some fundamental 
principles capable of providing reasonable and acceptable forms of modification, special quantities related to the 
cosmological or astrophysical observations are also expected. 

Particularly, the well known Noether's theorem is suitable for the purpose mentioned above, which has been 
extensively studied with allowed models given by the Noether symmetry approach, ie., in scalar field cosmology [26|, 
non-minimally coupled cosmology [27|, f(R) theory [28], f(T) theory [2{|, scalar-tensor theory [30| and quantum 
cosmology [31] , 

The complete solutions via Noether symmetry [32[ need to be confirmed, or expand it to General Noether symme- 
try [33[ . Physical contents of conserved quantities and their relations to cosmological or astrophysical observations are 
of great interest. Notice that not all of the conserved quantities derived from Noether's are physically meaningful, nor 
satisfy the astrophysical constraints, and there may also exist degeneracy phenomena between Noether symmetries. 

In this article, we investigate in detail several manifestations of the f(T) gravity according to Noether's theory and 
observational constraints. This piece of work is arranged as follows: In the next section we will briefly review the 
f(T) theory. The Noether symmetries for high-order f(T) manifestations are thoroughly researched in section three, 
ie., the exponential form and polynomial expression. Consistency of those models will be checked with astrophysical 
data-sets are also provided. In the last sections, we shall summarize some conclusions with discussions. 

II. BRIEF INTRODUCTORY TO /(T) 

In the teleparallel gravity, the vierbein field e^x^) plays the role of dynamical variables, which are defined as the or- 
thogonal basis of the tangent space at each point x^ in the manifold, namely, e^ • e^—^ , where r]ij=diag(l, —1,-1,-1) 
is the Minkowski metric. The vierbein vectors can be expanded in space-time coordinate basis: ej=e^<9 M , e*=e* dir^. 
According to the convention, Latin indices and Greek indices, both running from to 3, label the tangent space 
coordinates and the space-time coordinates respectively. The components of vierbein are related by e l e^=df and 

The metric tensor is related to the vierbein field ei(x fJ ') as 

Qiiv = Vijfal, (1) 

which can be equivalently expressed as rjij=g fJ-u e i ll e^ . The definition of the non- vanishing torsion tensor is given then 
by 

T% = e?(0„e* - d v e% (2) 

Evidently, T p vanishes in the Riemann geometry since the Levi-Civita connection is symmetric with respect to the 
two covariant indices. Differing from general relativity, teleparallel gravity relies on Weitzenbock connection, which 
is defined directly from the vierbein. 

In order to get the action of the teleparallel gravity, it is convenient to define the contorsion: 

^ = -~(T^-T^-T/-0. (3) 

Moreover, instead of the Ricci scalar R for the Lagrangian density in general relativity, the torsion scalar T describing 
the teleparallel Lagrangian density is defined by 

T = S^T"^, (4) 

where 

sr = \( K % + v tT ^ - v 1 ^)- ( 5 ) 

The modified teleparallel action of f(T) gravity is expressed as 

1 



16ttG 



■f(T), (6) 



where e=dct(e* )=y/— g- The teleparallel gravity can be obtained by setting f(T)—T. This modification is expected 
possibly to provide a natural way to explain the astrophysical observations, especially for the mysterious dark energy, 
as a motivation. 



III. NOETHER SYMMETRY FOR f(T) 

The field equations in modified gravity are partial differential equations (PDEs), and solutions of which can be 
provided and investigated by the Noether's theorem in a systematic and mathematical way, considered as the most 
elegant and systematic approach to compute conserved quantities, given by Noether in 1918. The conservation 
laws play a vital role in the study of physical phenomena so far. The Noether theorem states that any differentiable 
symmetry of the action of a physical system has a corresponding conservation law. The main feature of this theorem is 
that it may provide meaningful information regarding the conservation laws in the theory, ie., angular momentum can 
be explained by the rotational symmetry, and energy conservation due to the time translational invariance. Moreover, 
the Noether symmetries arc defined as transformations in the configuration space which preserve the Lagrangian 
form exactly. In the generalized symmetries, the configuration space is extended to the first or higher order of time 
derivatives of coordinates, which do not preserve the Lagrangian form exactly, except the modulo contact forms and 
exact differentials. 

Our present work is built on the base of flat FRW space-time, and the explicit description of torsion scalar T = —6H 2 . 
The method of Lagrange multipliers is adopted to set T as a constraint of the dynamics, we can obtain the point-like 
FRW Lagrangian 

L = a\f- f T T) - 6f T ad 2 . (7) 

The configuration space is {t,a,T}, so the Noether symmetry generator reads 

Q Q Q 

X = T (t,a,T)- + i,(t,a,T)— + <f>{t,a,T)—, (8) 

where r, tp and <f> are smooth functions of independent variables t and the canonical coordinates a and T. Then the 
first prolongation XM can be defined as 

lH = I + #,a,T)- + #,a,T)-. (9) 

The Lagrangian quantity should satisfy the following equation: 

X^L + (Dt)L = D B(t,a,T), (10) 

where B(t, a, T) is the gauge function, and D indicates the total derivative. The first integral is known as the conserved 
quantity associated with X by time integral, which is defined as 

fit fit 

I = B-TL-(ib-Ta)^=--(d-TT)^-r. (11) 

yY 'da v dT 

Expanding Eq. (fTOj) with the Lagrangian of f(T), a determined system of linear PDEs can be obtained: 

afrT a = 0, 

afrTT = 0, 

afripT = 0, 

1af T ^ a + f T ip + af TT (j) - af T T t = 0, 

-12af T i, t +a 3 (f-f T T)T a = B a , 

a 3 (f-fTT)r t = B T , 

3a 2 (f-f T T)iJj-a 3 fTTTcj) + a 3 (f-f T T)T t = B t . (12) 

In the following we will investigate several specific forms (models) of manifestation of f(T) theory with corresponding 
Noether symmetries, cosmological solutions and astrophysical data-sets constraints on relevant parameters of each 
model. 

A. Case I: /(T) ~ T n 
1. model depiction 

In this case where we assume that f(T) ~ T n , which is the simplest and commonly considered expression adopted 
in the f(T) theory [lj, [l8|]. According to the introduction above, taking the form to the determined system of linear 



PDEs (|T2|) . the solutions of r, </>, -0 an d B read 



C 1 t + C 2 , (13) 



= -~« 



iC 3 a- 3/2 ™ + Ci(l-2n) 



(14) 

-2Ci+C 3 a- 3/2n jr, (15) 

S = const, (16) 

where the index number n is assumed to be a real value, rather than an integer, and the Ci are the normalized 
orthogonal coefficients. By setting one of the Ci to 1 and the others to 0, we can get the corresponding Noether 
symmetry generators, which are deduced as 

d 2n - 1 d „ m d 
X 2 = |, (18) 

y 1 l-3/2n " , -3/2nrp " / in \ 

X 3 = --na _ + o T-, (19) 

where Eq. ()18[) apparently indicates energy conservation in the universe, while the physical meanings of Eqs. (|17j) 
and p9[) are a little more complicated and vague, which may suggest that time invariance being conjugated with 
cosmic inflation and the change of torsion. 

The commutator relations of the set {Xi t X2, X$} can be presented by the non-vanishing Lie bracket [Xi,Xj], 

2n — 1 

[X 1 ,X 2 ]=-X 2 , [X 1 ,X 3 } = X 3 , [X 2 ,X 3 }=0, (20) 

In 

thus the commutator relations of the set are closed. Moreover, the corresponding first integrals read 

h = a 3 T n t(2n-l)+4a 3 HT n - 1 n{2n-l), (21) 

I 2 = a 3 T n (2n-l), (22) 

h = -Aa 3 n 2 a- 3/2n HT n -\ (23) 

among which, I 2 and I 3 are time independent integrals. Setting ag = l (normalization of the scale factor at present), 
T = —6H 2 and 1 + z = - (relation between red-shift and the scale factor) conventionally, Eqs. (|2"2")l and (|2"3"|) finally 
indicate the same cosmological solution, which reads 

H = flo(l + z?' 2n , (24) 

where H$ represents the Hubble parameter at present. When n equals to one, an universe filled with non-relativistic 
matter only can be retrieved, according to J7 m = 1 for H(z) = Ho^Q m (l + z) 3 . From a much practical view, we can 
forecast that the observational data-sets constraints must favor n > 1 in order to provide explanations for dark energy 
and relativistic constituents (which include photons and neutrinos). 

The third solution (corresponding to symmetric generator X3) has been discussed in a previous research [291 ] . where 
the configuration space considered was incomplete. So the other two solutions (corresponding to symmetric generators 
X\ and X 2 ) are not found there. The I 2 can be regarded as the quantity of energy conservation according to the 
time translational symmetry X 2 = -Sr . The corresponding Hamiltonian E — 1 2 can also be derived by Eq. (l28l) , and 
the Eq. (|2"2"j) or (f2"3")l can also be analytically solved for a(i): 

8 ( t ) = (^)W3. t V3, (25) 

2n 

This solution is also obtained when Ii — 0, which means if we set the unknown conservation quantity I\ as zero, only 
one final solution (Eqs. (|24|) or ([251) ) can be found according to Noether 's theorem in this specific case. 

Through another direction, we can obtain the modified Friedman's equation by inserting the Lagrangian ([7]) into 
the Einstein's field equation, which read 

12H 2 f T + f = K 2 p m , (26) 

A8H 2 f TT H - f T {!2H 2 +4H)-f = p. (27) 



And the Hamiltonian for the from L = T(a, a, T, T) — V(a, T) of the Lagrangian 

■dL .dL T 

E = T^+a—-L 
dT da 

= T(a,a,T,f) + V(a,T) 

= -6f T aa 2 +a 3 (f-f T T) 

= a 3 (~l2H 2 f T ~f), 



(28) 



where T(a,d,T, T) = — Gfxaa 2 defines the kinematic energy term, and V(a,T) = a 3 (/ — JtT) denotes the potential 
energy term. Combining the first equation (|26p of modified Friedman's equation with Eq. (1281) , we obtain 



E = -a k p m => E = -k Q m pcr,o =*■ E = -3iJ o ri m , 



(29) 



where p m — Pmo/a 3 , £l m = Pmo/PcrO with p cr a — iH^/K 2 as the critical density and i?o is Hubble quantity at present. 
Combining the both equations above, we have 



3Hfo m = a 3 (UH 2 f T + f). 



(30) 



Then taking the form of f(T) —mT n T n for dimensional analysis, and a(0) = 1, the matter density parameter reads 



Q m = m (2n-l)a 3 (^) n , 
Jo 



which is useful to constrain the parameters n and m from the matter component tt r< 

1 L lri: 



(2n-l) 



(31) 



(32) 



2. Astrophysical data constraints 

Observational Hubble Parameter (OHD or H(z)) data-sets I34f , Baryon Acoustic Oscillation (BAO) data-sets 
(which contains data from 6dF 35], SDSS 36] and Wiggle Z 37] projects), and the latest Union2.f SN la (SN) 
data-sets [381 ] are taken for model constraining in our work, the best-fit results of parameters are listed in Table. |l] 
with corresponding minimal x-square value. Besides which the confidence ranges of parameter pair (n, h) are shown 
in Figs. Q] 
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FIG. 1: The lcr and 2a confidence ranges of parameter pair (n,h) (where h — Ho/(WO km ■ s _1 • Mpc -1 )), constrained by the 
observational H(z), BAO and SN la data-sets. The left contours show the confidence ranges constrained by H(z) data only, 
while constraints given by the combination of the three data-sets are shown as the contours on the right. 



TABLE I: Best-fit parameters of f(T) ~ T n model by H(z), SN and BAO with minimal x-square. 



data sets 


Kmin 


th 


(km/s/Mpc) 


n 


H(z) 


20.7488 




62.7074 


1.58407 


SN 


568.313 




69.1659 


2.34754 


SN+H(z) 


599.851 




68.8606 


2.20461 


SN+H(z)+BAO 603.965 




68.7674 


2.16819 



f{T) ~ T n 




FIG. 2: Contrasts between ACDM and T n model, where the shadow area is permitted by fitted parameters of various data-sets. 
For the bigger value of n, the H(z) curve becomes more deflected at high red-shift. Inconsistency is drown by the fact that 
smaller n fits data-sets better at high red-shift, while bigger n performs better at low red-shift. 

The best-fit results and confidence ranges provided by H(z) alone deviate significantly from the other data-sets, the 
best-fit parameter n and Hubble quantity h are smaller than those given by SN or BAO data-sets. The deviations are 
shown in Fig. [2] from which we find that this T n model predicts a static universe in the far future with H{z = —1) =0, 
besides it does not fit CMB observations as well, deviating greatly from the standard ACDM at high red-shift. It is 
very clear that /(T) ~ T n may not be an ideal approach of f{T) theory also because of its inconsistency as shown in 
Fig. [5] So other forms of expressions ie., f(T) ~ T m + T", should be taken into consideration. 



Case II: /(T) ~ ~£n T n and e" 



Since the previous model is not ideal according to astrophysical observation, we shall try other bold approaches 
which read f(T) ~ J2 n T n and e nT . Taking these forms to the determined system of linear PDEs (|T2|) . the solutions 
of r, (p, ip and B for both forms are respectively identical: 

(33) 
(34) 
(35) 
(36) 



(37) 



There exists only one corresponding Noether symmetry generator which reads 



T = 


Ci, 


il> = 


o, 


= 


o, 


B = 


const 
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x x 


d 
~ dt' 



The first integrals are derived as 



and 



7ia = 5Z[(2n-l)a 3 T"], (38) 



I lb = e nl a i (2nT-l), (39) 

where I\ a is the first integral of the symmetry of f(T) ~ J^ T", and ii b is the first integral of the symmetry of 

f(T) ~ e nT . The exact solution for /(T) ~ e nT is a ex: e' 2 ^-' 6 " and 7i & = 0, but H — const is apparently not suitable as 
a cosmological model which must fits observations. From the two cases discussed above we shall draw our conclusion 
that neither the simple T n model nor the extremely compact e nT model is acceptable theoretically and practically (by 
which we mean that any acceptable model should be both theoretically reasonable and observationally permissible). 
Since the two proposals have failed in deriving appropriate cosmological model discussed above, in the following, we 
will consider the truncated expressions of case II, like T" + T m or T + T n , which remain as another possibility. 

C. Case III: f(T) = aT + f3T n 

This case can also be regarded as a modification of case I by an addition of a linear term when considering the 
transition of the large scale of cosmological models to a smaller one, ie., the scale of solar system, for which we can 
also constraint the parameter n by corres pon ding conservation quantities. The form (T + T n ) has been researched 
with constraints from the solar system [39], |40(, and we will show its relation with Noether's theorem and cosmological 
observations. 

For f(T) — aT + (3T n , the only remained symmetry is the time translational symmetry and the corresponding first 
integral deduced from Eq. (|3"5|) as 



1 = aa 3 T + (3(2n - l)a 3 T". (40) 

at 

(2n - l)(-6H' z ) n - l (3H 2 + 2nH) + 'iH z + 2H = 0, (41) 



According to the conservation requirement §f =0, the Hubble quantity must satisfy the following condition 

-(2n-l)(-6F 2v ' ■'<•"'- ■>■■"' -" t2 
a 



which can be numerically solved as a cosmological model, we find the solutions are valid only when coefficient n is an 
integral number or zero, ie., n — ..., 2, 1, 0, — 1, — 2... (n=i is equivalent to n=l which can verified by Eq. (I41[) , which 

indicates that aT + ji\JT is equivalent to 7T with the same cosmological solution) , otherwise the value of the Hubble 
quantity given by the solution of Eq. (|4"Tj) is not real. 

Notice that, when the coefficient n = 0, the ACDM is retrieved. And we find out that f(T) = aT-\-(3T~ l is the best 
model of this case, which can explain and mimic the behavior of dark energy fairly well, and slightly deviates from 
the standard model as shown in the Figs. [3] for detail. The magnitude of the order of coefficient ratio j3/a is as high 
as 10 8 constrained by various astrophysical data-sets (the best- fit parameters are shown in Table. iHlwith confidence 
ranges given in Fig. [4]). 

When the magnitude of torsion T is small on large scales, the ft/T term plays a leading role in f(T), explaining the 
cosmological effect of dark energy. When the value of T is much higher on small scales, ie., in the solar system, the 
leading term is aT which is in accordance with the principle of general relativity. Besides, the constraining results 
predict a magnitude of the difference of T between its small-scale value and large-scale value, which should be no less 
than the order of 10 4 . 



TABLE II: Best-fit parameters of /(T) ~ aT + /3T" 1 model by H(z), SN and BAO. 
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FIG. 3: Contrasts of H — z relation between ACDM (the dashed line) and f(T) = aT + f3T~ 1 model with permitted parameter 
values (the shadow area). Reduction of cosmic age is inevitable since the Hubble quantity depicted by this specific f(T) 
modification is larger than the standard one as shown by figure on the right. 
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FIG. 4: la and 2 a confidence ranges constrained by various data-sets, where the solid lines represent constraints from SN 
alone, dashed and dotted lines represent constraints from BAO+H(z) and H(z) respectively. The overlapped region give a tight 
restriction for the parameter ratio /3/a between 1.7 x 10 and 1.9 x 10 . 



IV. CONCLUSIONS 



In the modified gravity theories, the explicit expressions of Lagrangian decide the exact solutions of the field 
equations and are important to gain insight into the physical content of the theories. While the Noether's theorem 
offers us a useful tool for researching the conserved quantities and symmetries of modified gravity theories, especially 
f{T) theory investigated in this article. 

As the above research, in the simplest case where f(T) ~ T", the full set of symmetries are expressed corresponding 
to generators X\, X2 and X3 as shown by Eqs. (|17 p -(|19 p . among which X2 is related to energy conservation or Hamilton 
condition. X\ is a special additional symmetry which contains tjL term. X3 is an additional symmetry, which may 
be related to the conformal symmetry. Nevertheless, the cosmological solution of symmetric generators X2 and X3 
possess the same form as Eq. ([M]). which can also be derived by setting Ii—O. Although the parameter n can be fitted 
by astrophysical data, the model f(T) ~ T n is not self-consistent and deviates greatly from the standard depiction 
of the universe. 



We also tried the exponential (e nT ) and polynomial (J2 n T n ) f° rms of f(T) theory, since the exponential (e nT ) are 
not theoretically acceptable (predicting a universe with constant Hubble quantity), and the expressions need to be 
truncated. 

Not only as a specific truncated form of case II, it is also meaningful to add T to the simplest form (f(T) ~ T n ) of 
case I to make it much more reasonable and suitable under scale transition, ie. from large scales which could be larger 
than 100 Mpc to the scale of the solar system. We specifically obtained aT + f3T~ l as best fitted model in our present 
work. And aT + f3yT is equivalent to 7T for the same solution. As shown above, astrophysical data-sets provide 
pleasing constraints on /3/a and the Hubble quantity as well, where the magnitude of the order of the coefficients' 
ratio ft /a. which has seldom been investigated in previous work can be also regarded as a restriction for the magnitude 
of torsion on different scales. When the torsion is small on large-scales, the 1/T term plays the leading role, explaining 
the phenomenon of dark energy. On the other hand, when T grows much large on small-scales, ie., in the solar system, 
the leading term is T itself which is equivalent to the general relativity. From the view of p henomenology, the critical 
magnitude between aT and f3T~ 1 is similar with the critical velocity of MOND theory [41|. The cosmic age predicted 
by this specific model is reduced by a tiny amount, by introducing higher order ^T" 1 term to the form of aT + f3T _1 
can alleviate that drawback. 

The time translational symmetry derived from Noether's theorem is mainly concerned in this work. Noether's theo- 
rem is also a feasible approach to discuss the conformal symmetry, which is used to discuss the Dark Energy by many 
modified gravity theories, ie., the Jordan frame and Einstein frame are connected by the conformal transformation. 

The methodology presented in this work can also be adopted for f(R) theory. Similarly, f(R) ~ R n can not satisfy 
the early- stage and the late-stage cosmic acceleration at the same time, which has been pointed by other related 
work [42J. Concerning a modified expression which resembles the work done in this article is inevitable. 
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